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Abstract 

We are presenting the basic ingredients of a stringy mechanism able to resolve both the 
Hagedorn instabilities of finite temperature superstrings as well as the initial singularity of 
the induced cosmology in arbitrary dimensions. These are shown to be generic in a large 
class of (4, 0) type II superstring vacua, where non-trivial "gravito-magnetic" fluxes lift 
the Hagedorn instabilities of the thermal ensemble and the temperature duality symmetry is 
restored. This symmetry implies a universal maximal critical temperature. In all such models 
there are three characteristic regimes, each with a distinct effective field theory description: 
Two dual asymptotically cold regimes associated with the light thermal momentum and light 
thermal winding states, and the intermediate regime where additional massless thermal states 
appear. The partition function exhibits a conical structure as a function of the thermal 
modulus, irrespectively of the spacetime dimension. Thanks to asymptotic right-moving 
supersymmetry, the genus- 1 partition function is well- approximated by that of massless 
thermal radiation in all of the three effective field theory regimes. The resulting time- 
evolution describes a bouncing cosmology connecting, via spacelike branes, a contracting 
thermal "winding" Universe to an expanding thermal "momentum" Universe, free of any 
essential curvature singularities. The string coupling remains perturbative throughout the 
cosmological evolution. Bouncing cosmologies are presented for both zero and negative 
spatial curvature. 
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1 Introduction 



Observational evidence strongly supports that during an early cosmological era, the matter 
content of the Universe was in (near) thermal equilibrium, with very high temperature. If the 
degrees of freedom are to be described by a set of local quantum fields, such a state results 
in a singular cosmology. Indeed, if we follow the cosmological evolution backward in time, 
using Einstein's gravity field equations, we are driven to the initial curvature singularity [l]. 
Even if a period of inflation preceded the high temperature phase, it is found in typical field 
theory models that the cosmological evolution begins at a singularity. 

In string theory we expect a drastically different picture to emerge since new purely 
stringy degrees of freedom can dominate the high curvature and high temperature regimes, 
leading to phenomena that do not admit a conventional field theory description [2], with 
Riemannian concepts breaking down. String oscillators and winding states become relevant 
around the Hagedorn temperature (which is of order the string scale M^), before the onset 
of curvature singularities, and drive a phase transition towards a new stringy thermal vacuum 
[3 12 . The simplest way to isolate the relevant critical phenomena is via the Euclidean 
description of the thermal system, where Euclidean time is compactified on a circle with 
period given by the inverse temperature ^}|7j- At temperatures just above Hagedorn, certain 
string states winding the Euclidean time circle become tachyonic. These instabilities can be 
lifted either by condensing the tachyons [5|[7], or by turning on special gravito-magnetic 
fluxes, which inject into the thermal vacuum non-trivial winding and momentum charges, 
as in 13 -16 . If a stable stringy phase gets realized, it could be that the back-reacted 



cosmological evolution is non-singular and the initial singularity is absent. 

A mechanism within which the Hagedorn instabilities of the string gas are resolved and 
the initial curvature singularity is bypassed was realized recently in a class of two-dimensional 



superstring cosmologies, the so called Hybrid cosmologies 16 . The scope of the present 
work is to show that the key ingredients of this mechanism are generic in a diverse class of 
higher dimensional superstring models as well. In all of these models, finite temperature is 



introduced along with non-trivial gravito-magnetic fluxes 13 -16 , which lead to a restoration 
of the thermal duality symmetry of the partition function: Z(/3//3c) = Z{l3c/(3). Here /3 
denotes the period of the Euclidean time cycle, attaining a critical value Pc at the self-dual 
point. At this critical point additional massless thermal states appear, enhancing the local 
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Euclidean gauge symmetry. Typical examples include the tachyon-free type II A/4 = (4, 0) 
models at finite temperature and in the presence of non-trivial gravito-magnetic fluxes, which 



are described in great detail in the literature 13 16 . The fundamental properties of these 
models, which can lead to the resolution of the Hagedorn and the initial singularity, are 



well understood from the recent study of the two dimensional Hybrid models 16 , and are 
exhibited below: 



• The canonical thermal ensemble is modified by turning on non-trivial gravito-magnetic 
fluxes, which lift the usual Hagedorn instabilities. The fluxes inject non-trivial winding 
and momentum charges into the thermal vacuum and render the mass of the would- 
be tachyonic states semi-positive definite. The tachyon-free models are equivalent to 
freely acting asymmetric orbifolds obtained by modding out with {—l)^^6o, Fl being 
the left-moving space-time fermion number and 5q an order-2 shift along the Euclidean 
time cycle. Essentially, the fluxes regulate the contribution to the free energy of the 
massive string states. 

• Not only is the resulting spectrum of thermal masses semi-positive definite, but also 
the partition function is duality invariant under /3 — )■ and it is finite for all values 
of /3. At the critical point, new massless thermal states appear extending the U{1)l 
gauge symmetry associated to the Euclidean time cycle to a non-Abelian [SU{2)L\k=2 
symmetry. This is a purely stringy phenomenon, absent in any conventional field 
theory model. The self-dual point /3 = /3c is realized at the so called fermionic point. 
This universal property of all such superstring models follows from the conformal 
transformation properties of the left-moving NS vacuum {hi = — \)- 

• The extra massless states at the critical (fermionic) point have non-trivial momentum 
and winding charges so that = ±1 and pr = 0. These two extra states together 
with the thermal radius modulus give rise to the SU (2) enhanced symmetry. At the 
critical point, the massless states give rise to non-trivial backgrounds which admit a 
localized brane interpretation in the Euclidean. 

• For /3//3c ^ 1, the asymptotic behavior of the thermal partition function is dominated 
by the light thermal momentum states giving rise to the characteristic behavior of 
massless thermal radiation in d dimensions, modulo exponentially suppressed contri- 
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butions from the massive string oscillator states: 

0{e-^/^^), (1.1) 



where n* counts the number of effectively massless degrees of freedom; is the Stefan- 
Boltzmann constant for radiation and Vd-i is the spatial volume. What is extremely 
important is that thanks to the thermal duality symmetry, the asymptotic behavior for 
/3//3c ^ 1 is dual-to-thermal, as it is dominated by the light thermal winding states: 

Here also, the oscillator states give exponentially suppressed contributions. The con- 
tribution of the massive oscillator states remains finite at the critical point, as the 
fluxes modify and effectively reduce the density of thermally excited massive oscillator 
states. In most cases, the contribution of the massive oscillator states never dominates 



over the thermally excited massless states due to asymptotic supersymmetry 13 - 18 



• This behavior indicates the appropriate, duality invariant definition of the temperature 
T, valid in both asymptotic thermal regimes. Defining the thermal modulus a by 
e°' = (3 /(3c, the duality invariant temperature is given by T = Tee"'""'. Thus the 
temperature in these configurations, and consequently the energy density and pressure, 
never exceed a critical value. The maximal critical temperature is given by Tc = 1/Pc- 
In both asymptotic regimes (T ^ Tc), the partition function can be expressed in terms 
of the self-dual temperature as follows: 

^ = rr (I) + O (e-^=/^) . (1.3) 

We conclude that the stringy thermal system has three characteristic regimes: The two 
dual phases of light thermal momenta and light thermal windings, and a third intermediate 
regime corresponding to the extended symmetry point, where vortices described by massless 
thermal states carrying non-trivial momentum and winding charges become relevant. 

The presence of the localized massless states is crucial since they can marginally in- 
duce transitions between purely momentum and purely winding states, thus driving a phase 
transition between the two asymptotic regimes. As in [l6j, this phase transition admits a 
geometrical description, in terms of a "T-fold", with branes localized at the critical point 
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gluing the "momentum" and "winding" spaces. This gluing mechanism was explicitly re- 



alized in the two dimensional Hybrid model 15 , 16 , where the partition function and its 
conical structure were determined beyond any a' approximation at the perturbative genus- 
1 level. In the Hybrid model, the ingredients described above not only treat successfully 
the Hagedorn transition, but also they lead to non-singular thermal cosmologies in contrast 
to field theoretic cases. In this work, utilizing the fundamental ingredients and especially 
the branes sourced by the extra massless thermal states, we show that non-singular string 
cosmologies also exist in higher dimensions. 

The plan of the paper is as follows. In section [2| we construct type II thermal vacua 
in arbitrary dimension c? > 2, which are free of Hagedorn instabilities due to the presence 
of certain gravito-magnetic fluxes. We show that up to exponentially suppressed terms, 
the corresponding partition functions are well-approximated by the contributions of the 
thermally excited massless states up to the critical point. In section [3} we present an effective 
action valid in all three regimes associated to the "thermal momentum phase" (/3 > /3c), the 
"thermal winding phase" (/3 < f3c) and the non-geometrical [SU{2)L]k=2 point (/3 = (3c). 
At the critical temperature, additional massless thermal states source negative pressure 
contributions to the effective action. Cosmological solutions free of initial singularities in 
arbitrary dimension are exhibited in section |4j They are compatible with perturbation 
theory throughout the evolution and describe bouncing cosmologies, where a phase transition 
connects a contracting thermal winding space-time to an expanding thermal momentum 
space-time. Solutions which are radiation or curvature dominated at both the very early 
and very late cosmological times are presented. Finally, our results and perspectives are 
summarized in section |5l 



2 Thermal duality and the Hagedorn transition 

In this section we construct tachyon-free thermal configurations, starting with type II A/4 = 
(4, 0) models in various dimensions. The left-moving worldsheet degrees of freedom give 
rise to 16 real supercharges, while the remaining right-moving supersymmetries are broken 
spontaneously via asymmetric geometrical fluxes [l3}jT6] . Geometrical fluxes generalize the 
Scherk-Schwarz mechanism |49j to string theory [6}|20||2T]. At certain points in moduli 
space, where the moduli participating in the breaking of the right-moving supersymmetries 
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attain values close to the string scale, the local gauge symmetry is enhanced to a non- 



Abelian one 13 16 . Finite temperature and quantum effects will stabilize these moduli at 



such extended symmetry points 22 24 . Interesting examples include the two-dimensional 
Hybrid models, where the right-moving sector is characterized by unbroken massive spectrum 



(boson/fermion) degeneracy symmetry (MSDS) [14 16 



Finite temperature is introduced along with non-trivial gravito-magnetic ffuxes, threading 
the Euclidean time cycle together with other cycles responsible for the breaking of the right- 



moving supersymmetries 13 -16 . These fluxes inject into the thermal vacuum non-trivial 
momentum and winding charges and lift the Hagedorn instabilities of the thermal ensemble. 
To see how this occurs, recall that for special values of the fluxes, the model is equivalent 
to a freely acting asymmetric orbifold of the form {—1)^'^6q, where Fl is the left-moving 



space-time fermion number and is a Z2-shift along the Euclidean time circle 13 . The 
genus- 1 partition function is given by 



Va-i f d\ 1 1 ^ , 5r p ra. 

(27r)'^-i 4^(^+i)/2 (r^r^)^ 2 ^ ^ > f)^ i (io-d,io-d)[5J 

2 a.h 



mo, no 

(2.1) 

where Vd-i stands for the volume of the large spatial directions. Here r(io_d,io-(i)[^] denotes 
the asymmetrically half-shifted lattice associated with the compact directions, leading to 
the breaking of the right-moving supersymmetries at zero temperature. The last line in the 
integrand gives the combined effect of finite temperature and the gravito-magnetic fluxes, 
amounting to a thermal action on the left-movers of the worldsheet. As a result, the orbifold 
deviates from the conventional thermal one in sectors with odd right-moving fermion number 
Fr. At the preferred right-moving extended symmetry points, these sectors are massive with 
masses of order the string scale. Sectors with even Fr and in particular the initially massless 
bosons and fermions are thermally excited as in the canonical ensemble. Thus the deformed 
model does not differ appreciably from the conventional thermal one at temperatures lower 
than Hagedorn. The presence of the gravito-magnetic fluxes however allows for the existence 
of stable phases above the Hagedorn temperature for which the canonical thermal ensemble 
fails to converge. 
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The left-moving Os sector carries non-trivial momentum and winding charges so that 



This is just enough to produce a state that is at least holomorphically or anti-holomorphically 
massless, ensuring the absence of tachyons in the spectrum of thermal masses for all values 
of the thermal modulus Rq. In addition, the model remains tachyon free under all marginal 
deformations of the transverse dynamical moduli associated with the compact manifold. The 
partition function is finite and invariant under thermal duality symmetry, Rq — > l/(2i?o) 
(together with the Ss ^ Cg interchange), with the self-dual critical point occurring at 
the fermionic point Rc = At Rq = Rc, additional thermal states become massless, 

enhancing the U{1)l gauge symmetry associated with the compact Euclidean time circle to 
a non-Abelian [SU{2)L]k=2 symmetry. Their masses are equal to 



At the fermionic point, the corresponding left- and right-moving momenta and the associated 
vertex operators are given by 



where Or are weight (0, 1) right-moving operators. Such states can marginally induce transi- 
tions between purely thermal winding and purely thermal momentum states, and in addition 
exchange the spinor chirahty, 5*8 o Cg |16|. As we will see, these states induce a univer- 
sal, non- analytic conical structure in the partition function Z as a function of the thermal 
modulus Ro, irrespectively of the dimensionality of the model. 

Thermal duality implies a maximal critical temperature and the existence of two dual 
asymptotically cold regimes dominated by the light thermal momenta and the light thermal 
windings respectively. The regime of light momenta, Rq ^ Rc, gives rise to a thermal phase 
with temperature T = where /3 = 27ri?o is the period of the Euclidean time circle. In the 
regime of light windings, Rq ^ Rc, the vortices can be interpreted by T-duality as ordinary 
thermal excitations associated with a large circle of period (3 = The corresponding 

temperature is given by T = = Thus the system at small radii is again effectively 

cold. The two phases are distinguished by the light thermally excited spinors: At large 
radii these transform under the S'g-Spinor of 50 (8), while at small radii they transform 



- min \p\ -pW = min (mo + -)(2no + 1) = 



(2.2) 




(2.3) 



(2.4) 
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under the conjugate Cs-Spinor. At Rq = Rc, we find the intermediate regime where extra 
massless thermal states appear and enhance the local gauge symmetry to a non-Abelian 



gauge symmetry. As in Ref. 16 , the extra massless thermal states give rise to genus-0 



backgrounds, present when Rq = Rc, admitting a Euclidean brane interpretation 16 . Since 
transitions between purely momentum and purely winding states can occur in the presence of 
such condensates, the branes "glue together" the light momentum and light winding regimes. 

Defining the thermal modulus a hj Rq/Rc = e'^ , the duality invariant expression for the 
temperature, valid in both the winding and momentum phases, is given by 

T^T,e-H r,^^-!-. (2.5) 

attaining a maximal critical value Tc at the self-dual point a = 0. As a result, the energy 
density and pressure in these models are bounded, never exceeding certain maximal values. 

To identify further universal features concerning the thermal effective potential, it is 
illuminating to compare two models of different dimensionality, where the mechanism leading 
to the resolution of the Hagedorn singularity is transparent: 



A d = 2 Hybrid model, where the first line of Eq. (2.1 ) is given by 



(-1)="' f r,«gl = ^ (f.. - s..) , (2.6) 



a, 6 



exhibiting holomorphic/anti-holomorphic factorization and right-moving MSDS structure. 



as exemplified by the identity V24 — 6*24 = 24. This model was analyzed extensively in 15 , 16 



• A (i- dimensional model, where the breaking of the right-moving supersymmetries occurs 
via the coupling to the right-moving space-time fermion number of the momentum and 



winding charges associated to a single factorized cycle, whose radius we denote by Rg 13 

1 1 f ,^d+b ^[5]"^ p ran _ ^(9-d,9-d) 



a, 6 



mg,ng 

(2.7) 

The breaking of the right-moving supersymmetries has a similar algebraic structure to the 



temperature breaking that acts on the left-moving characters, see Eq. (2.1). The Euclidean 



partition function is invariant under the Rq -H- Rq exchange 13 . The fermionic point 
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Rg = Rc corresponds to an extended symmetry point, with the U{1)r gauge symmetry 
associated with the X^-cycle getting enhanced to [5'f/(2)^]fc=2- As we aheady stated, finite 
temperature and quantum effects give rise to an effective potential that stabihzes dynamically 



the value of the Rq modulus at the extended symmetry (fermionic) point 22 , 23 . This 
behavior is drastically different than that of the thermal Rq modulus, where the back-reaction 
in the Lorentzian drives a cosmological evolution towards smaller temperatures. All other 
spectator moduli are either stabilized at extended symmetry points or frozen at values of 



order unity 23 , 24 . The string coupling is taken to be sufficiently weak. As we will show, it 
remains smaller than a critical value during the induced cosmological evolution. Despite the 
lack of MSDS structure in the higher dimensional cases, the massive states are characterized 
by asymptotic supersymmetry, thanks to the asymmetric nature of the left- and right-moving 
supersymmetry breakings. This fact also explains the absence of tachyons from the spectrum 



of thermal masses 13 17 



In the two dimensional Hybrid model, the genus-1 partition function takes a very simple 
form, which makes transparent the conical structure of the thermal partition function. The 
expression can be derived beyond any a' approximation, thanks to the unbroken MSDS 



symmetry characterizing the right-moving sector, and it is given by 15| 16 
'Z'Hybrid - . / ~ 1 



Vi V 2Ro 



24 X i?o + 777^ - 24 X 



Ro 



2i?n 



24V2e~l'"l. (2. 



The essential feature is a discontinuity in the first derivative of ^Hybrid as a function of the 
thermal modulus a, signaling a phase transition between the two dual "momentum" and 
"winding" regimes. The discontinuity occurs at the self-dual, fermionic point a = 0, and it 
is sourced by the 24 complex lowest mass states in the 08V24-sector, which become massless 
precisely at this point. Due to the unbroken MSDS right-moving structure, there are exact 
cancellations between the massive fermionic and massive bosonic oscillator states, and so 
both regimes at \a\ > comprise phases where the equation of state is effectively that 
of massless thermal radiation in two dimensions. Despite the cancellations in the massive 
sector, stringy behavior survives, with states carrying both momentum and winding charges 
becoming massless, inducing a phase transition at the critical point a = 0. 

In terms of the duality invariant temperature, the partition function is simply given by 

%^= 487rT, (2.9) 
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attaining a maximal value at the critical temperature T^. The existence of this maximal 
temperature is the crucial difference from field theory thermal models where the temperature 



is unbounded. The numerical coefficient in Eq. (2.9) is easily understood as follows. In 
d dimensions, the partition function corresponding to massless thermal radiation can be 
written as 

-^ = n*S, r^-\ (2.10) 

Vd-l 

where n* is given in terms of the numbers of the initially massless bosons and fermions and 
Erf is the Stefan-Boltzmann constant: 

«-=""+«"'^. 2.=^i:^cM. (2.11) 

Taking into account that initially there are 8 x 24 massless bosons and 8 x 24 massless 
fermions in the two dimensional Hybrid model and that S2 = 7r/6, we obtain n*^Ti2 = AStc. 

We now proceed to analyze the (i-dimensional model. The non-analytic, conical structure 
of the partition function occurs at the fermionic point Rq = Re- We can identify it by 
computing Z for i?o > Rc and also for Rq < Rc, and then utilize thermal duality to connect 



the two regimes so as to obtain an expression valid for all radii (see Refs 15 , 16 for detailed 
discussions concerning the Hybrid model and also [9j for the two-dimensional Heterotic 
strings). For Rq > Rc, we Poisson resum over the momentum quantum number ttiq, and 
map the integral over the fundamental domain to an integral over the strip |25| involving 
the (mo, no) = {2k + 1, 0) orbits only. In the strip representation, the winding contributions 
and in particular those of the Og and Cg sectors to the fundamental domain integral, are 
mapped to Vg- and S'g-sector momentum contributions in the UV region of the strip, T2 < 1. 



Mapping the integral over the fundamental domain to an integral over the strip as in 15 , 16 
(see also [25]), we obtain 



00 „ 

1~ J\\ 



Vd-i ^t;;yii4(2vr)^-vr/^ V'^ 



k=0 



R9=Rc 

(2.12) 



where the moduli associated with the T(^g^d,9-d) lattice are taken to be of order unity. 

For large Rq, we split the contributions to the integral in two pieces: (i) the contribution 
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of the thermally excited massless bosons and fermions and (ii) the contributions of the 
massive states. 

(i) Massless contributions: 
The contribution of the initially massless bosons and fermions is given by 



= ^R>y:r ^^^W e-<^'«<'^ 16 X (8 + 2), 



' massless 

fe=0 



- 8 X (8 + —— - n l.d- 



2d-l pd-1 p'^-^ ' 

(ii) Massive contributions: 
To determine the contribution of the massive states, we must first compute the integral over 
Ti, which imposes level matching. The integral over T2 gives the final result in terms of 
Bessel functions. Due to the presence of the gravito-magnetic fluxes, there are alternating 
signs between the contributions of states in the Fr-cvcu and F^-odd sectors (right-moving 
bosons and right-moving fermions). The massive contributions become 

where the first sum is over individual degenerate boson/fermion pairs, having the same right- 
moving fermion number. Thanks to the (—1)^^ alternating signs, the effective density of 
states gets reduced drastically, as compared to the canonical thermal ensemble. This is a 
signal of right-moving asymptotic supersymmetry, replacing the exact right-moving MSDS 
structure of the Hybrid models. The massless sector on the other hand contributes identically 
as in the canonical thermal ensemble. For (3 ^ /3c, the arguments of the Bessel functions are 
large, leading to an exponential suppression of the contributions of all massive states. In this 
regime, we recover the characteristic behavior of massless thermal radiation in d dimensions: 



' massless 



By thermal duality, the behavior for f3 <^ is dual-to-thermal, yielding 

d-l 



(2.15) 



^ ^^^'^ ^ ^ > (2.16) 



This T-dual result for Rq <^ Rc can also be obtained if we first Poisson resum over the winding 
number no and utilize the momentum-unfolding to map the integral over the fundamental 
domain to an integral over the strip. 
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To complete the discussion, we have to analyze the behavior of the partition function in 
the intermediate regime when Rq is close but still larger than the fermionic point Rc. In 



this case, we must examine the contribution of the r2 — )■ region to the integral (2.12), 
due to the exponential growth in the density of massive states in each sector of definite 
parity separately. The individual contributions of the massive states are larger in this 
region. However, as we will see the contribution from this region is drastically reduced due to 
right-moving asymptotic supersymmetry, whose origin is the insertion of the (—1)^^ phase. 
This property also explains the absence of physical tachyons from the spectrum of thermal 
masses. The asymptotic r2 — )■ oo region is dominated by the lightest string states giving 
rise to the thermal massless radiation contribution. As we show below this contribution will 
turn out to be the dominant one. 



To proceed further we need to determine the r — )■ limit of the integrand in Eq. (2.12). 
To this end, it is convenient to rewrite the integrand in terms of shifted lattices 



oo „ 



e ^2 



where 



x^^^(r[119l-r[;i9l-r|;]9l)U.^^. 

r|](B) = « ^e-'#l«*+-'>+'="+''W'. (2.18) 



m,n 



We then apply the modular transformation t ^ t = —l/r to the following expression 



1 



~d/2-l 



6i ^{9-d,9-d) 1 

^4 {rif]Y fj 



(2.19) 



which appears in the integrand of Eq. (2.17). The last expression can be expanded in powers 
of g = e^^'^, in the limit g — 0. For Rg = Rc, we obtain the following leading terms 

- 8 ) f 8 + 2 + 2 g3 ^^'A , (2.20) 



~d/2-l 
'2 



Essentially, only the left /right level matched terms contribute due to the integration over ri. 
Keeping these only, we get 

^ : (8 + 2) = T^^'^'^S X (8 + 2) . (2.21) 



~d/2-l 
'2 
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The absence of exponential growth in this factor is the signal of asymptotic supersymmetry, 
as was already stated before. Notice that for d > 2, this factor goes to zero as a power law. 
So, the contribution from this region is estimated to be 

where t* = ^ is a sufficiently large cutoff, L ^ 1. The overall contribution is thus exponen- 
tially suppressed. 

Therefore the behavior of Z for i?o > Rc is controlled by the thermally excited initially 
massless states everywhere and up to the critical point. The same conclusion can be reached 
in the regime Rq < Rc by thermal duality. Gluing the two regimes in a duality invariant 
way gives 



modulo the exponentially suppressed contributions in the three effective field theory regimes. 
This result also implies that in each of the two thermal phases the various thermodynamical 
quantities enjoy the standard monotonicity properties as functions of the temperature, with 
the specific heat being positive up to the critical point. 

To conclude, the above result reveals a universal conical structure at cr = 0, irrespectively 
of the dimensionality of the model. All of the above manipulations, including thermal duality. 



amount to approximating the last line of equation (2.12) with a factor of order unity: 

2i?,el"l V / —jj^ exp ^ ^ e^l"! 8 + 2 . 2.24 



Thanks to the analytic properties of the left-moving characters, only the massless level 



contributes in the above integral. This is the generalization of the Hybrid model result 15 16 
to arbitrary dimensions via right-moving asymptotic supersymmetry. 



3 Effective action (s) up to genus- 1 



From the thermal effective field theory point of view, there are at least three different effective 
actions associated to three possible a'-like expansions, each being valid in one of the three 
characteristic regimes. Namely: 



Rq 
Rr 



>1 (a>0). 



Rc 
Ro 



>1 (-a>0). 



Rq 

Rr 



Rc 
Rn 



< 1 (a ~ 0). 



(3.1) 
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Taking into account the behavior of the thermal partition function in each regime, the 
thermal duality symmetry, as well as the branes which glue the dual "momentum" and 
"winding" phases, we will construct a d-dimensional effective cosmological action valid in all 
regimes simultaneously and derive the associated equations of motion. In the Lorentzian, 
the branes are spacelike, appearing at any time when the temperature reaches its critical 
maximal value T^.. As we will argue, these source localized negative pressure contributions 
to the effective action. 

In the two asymptotic regimes (cr — >■ ±00) dominated by the light thermal momenta and 
the light thermal windings respectively, the effective action admits the well known sigma- 
model descriptions, defined via the corresponding a'-expansions. Thanks to thermal duality, 
both asymptotic regimes can be simultaneously described by a unique expansion in terms 
of the duality invariant temperature T — Tce~^'^^. All thermodynamical quantities such as 
the temperature, the energy density and the pressure are given in terms of manifestly du- 
ality invariant expressions involving the absolute value of the thermal modulus \a\. In the 
Euclidean, the regime \a\ — )■ is well described in terms of the [SU{2)L]k=2 CFT associated 
with the fermionic extended symmetry point. In this regime, we have to include the con- 
tributions of the extra massless thermal states, responsible for the phase transition, both at 
the genus-0 and genus-1 levels. The genus-0 contributions admit a brane interpretation with 
a tension determined by the allowable non-trivial backgrounds of the extra massless thermal 
scalars, df^ip^ ^ 0, where the gradients are along the directions transverse to Euchdean time. 

During the cosmological evolution, the thermal modulus o acquires non-trivial time- 
dependence, o"(t). Since all fields are functions of |cr(r)|, their second time-derivatives may 
give rise to localized singular terms proportional to 



where the temperature reaches its critical value at times (i = 1, . . . , n) so that o(t:^ = 0. 
As we will see, these singularities are naturally resolved by the presence of the spacelike 
branes, localized at the times r^. 

The relevant representations of the winding-like field theory are space-time left-moving 
Vectors Vg and space-time anti-Spinors Cg- On the other hand, in the momentum-like field 
theory, the relevant operators are the left- moving Vectors Vg and space-time Spinors Sg. 




(3.2) 
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At the branes, the theory is self-dual and both the Spinor and anti-Spinor representations 
coexist together with extra massless states with non-trivial momentum and winding charges, 
{Pl,Pb) = (il) 0), triggering the transition of the winding-like field theory based on Vg — Cg 
to the momentum-like one based on Vg — Sg. 

The above ingredients lead to an effective c?- dimensional dilaton-gravity action (up to the 
genus-1 level), able to describe simultaneously and in a consistent way the three regimes: 

5 = iSo + iSi + iSbrane, (3-3) 

where 



2 

Si= I d^x^P, (3.4) 



'5brane = -^ J d'^X^/g^C '^'^ Ki 5 {t - T,i) 



- Sq is the genus-0 dilaton-gravity action written in the string frame. 

- Si is the genus-1 contribution of the thermal effective potential —P. 

- '^brane IS the spacelike brane contribution at the phase transition giving rise to localized 
negative pressure. It is sourced by the additional massless scalars ip^ from the OgVg and OgOg 
sectors at the extended symmetry point a = 0. These extra massless states parametrize a 
manifold, which up to discrete identifications, takes the form of a coset space, 

of real dimension 2q which depends on the number of massless states coming from the right- 
moving (non-supersymmetric) sector. In a suitable parametrization, the associated Kahler 
potential is given by 



= - In (^(r° + F°)^ - (F° + F")^ j , a = l,...,g-l, (3.6) 

with the metric given in terms of the holomorphic and anti-holomorphic derivatives of K: 

ds^M^") = dad^K dV'dY^ = K^-p dV^dY^ . 

The smallest possible value for g is 10 — d, which occurs when there is no extended right- 
moving gauge symmetry. In the models of 13 described in detail in section [2| q = 2+(l0—d), 
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with the gauge symmetry of the right-moving sector being extended to I-Lr = SU (2) x U{1)^~'^ 
of dimension 3 + {9 — d) = q. The maximum value for q is 3(10 — d), occurring when the 
gauge symmetry of the right-moving sector is extended to Hr of dimension 3{10 — d). In 
the simplest models with this property, Hr = SU{2y^~'^. 



3.1 The S-brane action 



The microscopic origin of the brane term in Eq. (3.3) follows from the underlying description 
of the system at the extended symmetry point. We will be interested in obtaining homoge- 
neous and isotropic solutions of the bulk. The scalars give rise to a tree-level localized 
action 

5brane = -j dcT d^^^ X ^ c'^'*' g''' G u f,^' d,^' 5{a) , (3.7) 

where /i = l,...,(i— l,(yf_L = det gf^y and Gjj (or K^^) is the metric in the field configuration 
space A^^''. The equations of motion of the scalars i-p^ take the form 

29^(e-2^v^ g'^'' Gijdo^') - e'^-^v^ g'^\diGKj)df,v''d,^' = . (3.8) 



Our aim is to establish that these equations admit non-trivial solutions which are consistent 
with the homogeneity and isotropy requirements and yield the localized brane contributions 



in the effective action as described by Eqs (3.3) and (3.4). To this end, it suffices that at 



each instant {i = 1, . . . ,n), when the temperature reaches its critical value, the induced 
metric be proportional to the spatial metric: 

h/iu = Gij df,^^ dc^-^ = 9[iu, (3.9) 

where the positive constants. When this happens, the stress tensor of the scalars is 

consistent with the symmetries of the spatial metric, and therefore with homogeneity and 



isotropy. Moreover, the action (3.7) takes the familiar form of the Nambu-Goto action for 
branes. 



5brane = - ^ Ki d X 6 .Jgl 5(r - Tj) 

i 

= -5^(rf-l)^«:J^ y d'^xe-'"^ ^det {Gijdf,4>'do4>J) 8{t - r,) , (3.10) 



where we have used Eq. (3.2). Thus, Kj is interpreted as a brane tension. 
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To exhibit the solutions, we first consider the following embedding of space fl'^ ^{x^) into 
the field configuration space, n°'~^(a;^) — )■ M.'^^{(p^) : 

dp_(p^ = 5^, /t, / = 1, . . . , c? — 1, and = const., / = d, . . . , 2g. (3-11) 

The embedding exists provided that the dimensionality of the scalar field manifold is bigger 
or equal to the dimension of space: 2g > — 1. It also implies that the induced metric 



satisfies: /i^p = Gj^^. Thus Eq. (3.9) imposes that the spatial metric gj:^c and G^p are 
isomorphic 

Gfiu = ^ _ ^ dfiu ! (3.12) 
upon the identification Lp^ = x'^. Under the above circumstances, the field equations of 



motion (3.8) for I = fi < d — I become: 

29^(e-2'^v^) - e-2<^v^ g^'df.goa = =^ 9/^0 = 0, (3.13) 
consistently with the homogeneity of the dilaton field (p, and where we have used the identity 

g"^df,gya = 2df,\ogy/gI. 



Eqs (3.8) for I > d — 1 must also be satisfied. The above discussion makes it clear that 
the geometrical structure of the field manifold A4'^'^{(f^) is crucial, since it constrains the 
possible embeddings of Q'^~^{x^) into Ai'^'^{Lp'). We are mainly interested for the isotropic 
embeddings of the hyperbolic space H'^^^ and the flat space F'^~^ , with curvature and metric 
given by 

(^) H''-^ ■.k = -l (for > 2) , g(,, = 5^,, (ii) F''-^ : A; = , g^^o = ainf , 

\x ) 

(3.14) 

where a(Tj) is the scale factor at time Xj. 

{i) Hyperbolic embedding i^*^"^ — > M."^^ 

The hyperbolic embedding turns out to be naturally realized thanks to the geometrical 
structure of the field manifold Ai'^'^. Indeed it is sufficient to utilize the sub-manifold 
Ki{y^) C M^'^iY^) with = Re (J = 0, . . . , g - 1) non trivial and uj^ = ImY^ fixed. 
The sub-manifold K'^{y') naturally contains the desired if^^ factor: 

K'^ ^ H^-' X 50(1, 1) = X soil, 1) . (3.15) 
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The metric on K'' follows from the Kahler metric K^^ and takes the form: 



2^2 ^2 



ds\K'^) = + ^ with r^ = (yy-(y^y. (3.16) 



The imaginary parts are frozen consistently with all equations of motion. The metric of 
the (g — l)-dimensional hyperboloid H'^~^ is obtained by writing 



-{dyy + {dy'Y = -dr' + r' {dH'^^y , (3.17) 

which shows the explicit factorization of K'^: 

ds'^iK'^) = ^ {dC + {dH'^-^f) , C = Inr . (3.18) 

Fixing further the field C to be constant, as allowed by the equations of motion, the em- 
bedding of the spatial hyperboloid H'^~^{x'^) into H'^~^{ip^) is automatic. This can be done, 
provided that q> d, ip^ = and the extra q — d fields ip^,i = d, . . . q — 1 are also frozen as 
allowed by the equations of motion: 

d^ )- (y,d-l)2 + (y,d-l)2 =^ ds {H )- (3.19) 

When q = 2 + (10 — d) as in the models of |13|, the constraint d < q implies d < 6, while 
when q takes the maximal value 3(10 — d), we must have d < 7. These embeddings satisfy 



the relation (3.9) with the tension /tj fixed in terms of the scale factor: Hi = {d — l)/2a(rj)^. 
Such a tension leads to highly curved cosmological solutions with Ricci curvature of order 
one, as the value of a(rj) is fixed to be of order one by the cosmological equations (see section 



4.2j). 

To relax the constraint on the tension Kj, we take into account discrete identifications 
in the field configuration space in order to obtain solutions with a non-trivial wrapping 
number. With suitable discrete identifications, the submanifold parametrized by the fields 
(f^{I = l,..., d—1) becomes a finite volume hyperbolic space of the form H'^~^/r, where F is 
a subgroup of the discrete duality group S0{2, q; Z) of Ai'^'^. The two dimensional examples 
correspond to the familiar higher genus Riemann surfaces. Homogeneous compact hyperbolic 
manifolds in dimension > 3 are characterized by the property of rigidity, which implies that 



there are no massless shape moduli 26 , and are locally isotropic. The volume is determined 
by the radius of curvature L and the topology of the manifold: Vol(if°'~^/F) = L°'~^e°, where 
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a is a constant determined by the topology and L ~ 1 for the metric on the field subspace. 
The topological factor e" is unbounded from above. Taking the spatial manifold Q'^~^{x^) 
to also be a compact hyperbolic manifold of large volume (and suitable topology allows 
for embeddings with arbitrary wrapping number w. Consequently the brane action is finite, 
given by 

^brane = " e'^-^^ Vol(fi^-l) = -K, [ dx'~h~''^^ ^ , (3.20) 

implying that the brane tension is = c* w'^/a{Ti)'^. Here c* is a factor determined by the 
topology. The wrapping number w being arbitrary, the tension can be kept arbitrary. 

(ii) Flat embedding F'^^^ — > M'^'^ 

A way to realize the fiat embedding is by utilizing a (rf — l)-dimensional flat section F'^~^{ip^) 
of Ai'^'^. The isotropic embedding of fiat space F'^^^{x'^) into the fiat section F'^~^{ip'^) is 
defined by 

2/?i 
d 



du^l- = \l ^ a{n) 5^ , (3.21) 



giving rise to Eq. (3.9) and the brane action with tension Kj. We display below examples of 
such fiat sections. 

• For = 2, we can utilize the S0{1, 1)^ factor of the submanifold K'^ parametrized by the 



field C in Eq- (3.18) to carry out the embedding — )■ A^^'^. All other fields are frozen 
consistently with the equations of motion. 

• For c? = 3, we utilize the 5*0(1, l)(j x 50(1, 1)^ submanifold of K'^ parametrized by the 



fields C and ^ = ln(^"'-\ see Eq. (3.19). 



• For > 4 we utilize the fiat section of H'^^^{u^) obtained in the limit u'^^^ — )■ oo, together 
with the 5*0(1, 1),^ factor. To this end, we define the rescaled fields by setting = u'^~^(p^ 
(i = 1, . . . , d — 2), in order to obtain: 

I [de + idH'-r] = I [dc + id^r] + J^yd^^du^-^ + ^^^)^idu'-'r 

= \ [de + [d^r] + o{^) . (3.22) 



The field C, and the d — 2 fields v?* are utilized to realize the fiat embedding of Eq. (3.21). 



The 2q — {d — 1) extra fields are frozen, consistently with the equations of motion (3.8) 



^See e.g. [27j[28j for discussions concerning this possibility in the context of cosmological and other 
appHcations. 
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including that of u'^ ^ in the hmit u'^ ^ — )■ oo. The reahzation of this isotropic embedding 



imposes the constraint d < q. When q = 2 + (10 — d) as in the models of 13 , this implies 
d < 6, while when q takes the maximal value 3(10 — d), d < 7 is required. 

Finally let us note that isotropic embeddings are also possible when space is isomorphic 
to the ((i — l)-dimensional sphere S'^^^. These embeddings, the restrictions on the dimension- 
ality of spacetime and the structure of the bulk solutions are currently under investigation. 
This completes our discussion about the origin of the brane contributions in the effective 
action. 



3.2 Equations of motion 

Looking for homogeneous and isotropic cosmological solutions in dimensions d > 2, the 
dilaton field is a function of time only and the metric 

ds^ = -N^rfdr^ + a(r)2(rffi'^-i)2, = ^(r) (3.23) 

involves the line element {dVf^'^Y of the {d — l)-dimensional Einstein space with curvature 
k. To derive the equations of motion, we utilize the analytic expression for the Ricci scalar 
curvature 



^_ 2(^-1) 
7V2 



■a id -2) f^2.kN^\ rrN 



where H =-. (3.24) 
a 



The pressure P is determined by the genus- 1 Euclidean path integral Z (the thermal partition 
function) as 

P = T,e-H^^. (3.25) 
Vd-i 

It is important to note that in the sigma-model frame, P is a function of the thermal modulus 
|cr| only and there is no dependence on the dilaton field 0. 

The above considerations lead us to the following equations of motion: 

(i) N -equation 

]^{d-l){d- 2) (^H^ + k^^ = 2{d - 1)H<P - 202 + e^'^N^p, (3.26) 
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where the energy density is given 



P^-P-l^^-r.e-HJLM). (3.27) 

d\a\ d\a\ \ Vd-i J 

(ii) a-equation 



{d-2)^ + i(d - 2)(d - 3) (h' + k^] -{d- 2)H^ 



+ 2{d - 2)H^ - - 20^ - e^^N^P + ^ KiN 5{t - n). 



(3.28) 



N 

(iii) (f)- equation: 



0^- 










d-2 


\ a 


2 



(3.29) 



T - Ti 



It is useful to disentangle the second derivatives of the scale factor and dilaton field in the 
last two equations. Of particular interest is the linear combination which leads 

(ii)' trace equation 



2 \ a 



^ ^^"-if|j-^(rf-l)(d-2)(/72 + A:^)=e2^iV2[p-(rf-l)P] (3.30) 



and shows that the second derivative of the scale factor is always finite. In particular, this 
implies that a is continuous, even at the branes localized at r = Xj. On the contrary, another 
equation for the dilaton is 

(iii)' (p-equation modulo trace equation 
20-402 + 2(rf-l)i70-20^ = e2'^Ar2[(rf-l)P-p] -^^,^N 5{t - n), (3.31) 

i 

which shows that the first time- derivative of the dilaton is discontinuous across r^. This 
discontinuity is resolved by the presence of a spacelike brane, whose tension must satisfy 

.. = 2fa^. (3.32) 



2To derive p, we make use of ^'^^^'i^gy^'^'" - 5{r' r)P + A^|^^ ^here ^ = as 
follows from the fact that /S^el'^^^)!^ = N{t). 
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This resolution of the discontinuity via branes provides a novel mechanism in obtaining non- 
singular bouncing cosmologies, which in fact remain in a perturbative regime throughout the 
evolution. Such dynamical behavior is induced by the presence of spacelike branes providing 
localized negative contributions to the pressure 

PB = -^e-2<^(^")/€,5(r,), (3.33) 

i 

thus evading the conditions set by the standard singularity theorems on realizing non-singular 
bouncing cosmologies (l]. 

Before solving the above cosmological equations, it is important to stress that these 
give rise to an integrable relation of fundamental interest, namely the Entropy conservation 
equation: 

(s) Entropy equation 

{p + P) + {{d - l)H + \a\) {p + P) = 0. (3.34) 

Integrating once, we obtain a quantity 5*, which we may interpret as the conserved thermal 
entropy in a comoving cell of volume a'^~^: 

(s)' Entropy conservation 

"j^iP + P) = S. (3.35) 

Since the energy density p and the pressure P are bounded, attaining their maximal values at 
the critical temperature Tc, the scale factor a is bounded from below, acquiring its minimal 
value at the critical point. Thus the big-bang singularity a = of general relativity is avoided 
in all such Hagedorn-free string models. The constant thermal entropy S can be computed in 
the asymptotic regime where the system is radiation dominated and the dilaton is constant: 



p ~ (d- 1)P ~ (rf- l)n*SdT'^ as |a|^oo =^ S = d-ft^^ n*^d. (3.36) 

where the constant 7oo is given in terms of the asymptotic values of a and T, = 
lim|^l_,oo oT. 

As we demonstrated, the stringy thermal system is effectively radiation dominated in all 
regimes. As a result, an important consequence of the applicability of the entropy conser- 
vation in all regimes is that 7 = aT is almost constant everywhere and up to the critical 
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temperature. The largeness of the entropy observed at late cosmological times implies that 
the size of the Universe at the critical point is already large. This fact guarantees the validity 
of our perturbative approach as we will see later, and also, the connection of the so called 
entropy and oldness problems of standard Big/Bang cosmology. 



4 Stringy non-singular cosmologies 

In this section, we exhibit non-singular cosmological solutions in various dimensions. We will 
show that the mechanism for the resolutions of both the Hagedorn and the initial singularity 
problems are generic in all space-time dimensions. This result follows from the universal 
properties of the partition function discussed in detail in sections 2 and 3. Namely: 

• The thermal partition function has a conical singularity as a function of the thermal 
modulus a at the critical point cr = 0, irrespectively of the space-time dimension. This 
implies a phase transition between the light-momentum and light winding-effective 
field theory regimes. At the phase transition, spacelike branes appear connecting the 
two asymptotic regimes. 

• Across the branes, the behavior of the scale factor and dilaton field is such that a is 
continuous, while there is a discontinuity in related to the brane tension, see Eq. 



(3.32). 



Two-dimensional examples have already been exhibited in Ref. 16 in the framework of the 
Hybrid models, where the contributions of the massive modes to the partition function Z 
cancel exactly thanks to the MSDS structure characterizing the right-moving sector. In 
the higher dimensional models presented in section |2} the exact structure of the partition 
function is more involved. However, knowing the asymptotic behavior of Z as cr — t- ±oo, as 
well as the asymptotic density of the right-moving states, allow us to determine the essential 
thermodynamical properties for all temperatures T < Tc. Indeed, thanks to the asymptotic 
right-moving supersymmetry characterizing the models, we can well-approximate Z with 
the contribution of the thermally excited massless states in both the a > and the a < 
regimes, as was explained in section [2j Thus, in what follows we set 

p = {d-l)P= {d - I) rCT.dT'^ = (c/- l)n*Sdr^'^e-'^l"l. (4.1) 
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This last expression captures the essential conical structure of the higher dimensional cases, 
generalizing the two dimensional Hybrid results. It allows us to obtain analytic cosmological 
solutions via spacelike branes, as in the two-dimensional case. The conical structure arises 
from matching the dual momentum and winding regimes, which admit distinct effective field 
theory descriptions. The intermediate a = regime gives rise to the spacelike branes. 

The brane interpretation gives us the possibility to look for branches of solutions defined 
in time intervals (rj,rj+i), with o"(rj) = 0. These intervals are connected to each other by 
spacelike branes localized at Tj. From the world-sheet point of view the transition occurs 
via a condensation of massless thermal states carrying non-trivial momentum and winding 



charges associated to the vertex operators 0± defined in Eq. (2.4). The existence of these 
additional marginal operators at cr = gives rise to the transition between the thermal 
winding states at cr = 0_ and the thermal momentum states at cr = 0+ and vice-versa. The 
Universe may experience a series of such phase transitions. For a cosmological evolution to 
be consistent, the following constraints must be fulfilled: 

1. The dilaton and scale factor must be continuous across the branes. 

2. The scale factor must be bounded from below, a(r) > Oc, since T oc 1/a cannot reach 
values larger than the maximal temperature T^. 

3. Across any brane, 4>{Ti-) > 0(Tj+) is required, since the discontinuities in the first time 



derivative of the dilaton are resolved by a positive brane tension (see Eq. (3.32)). 



4. Since the first derivative of the scale factor is smooth at the branes, with the scale 
factor reaching its minimum value, we must have d(rj) = 0. 

5. To maintain a perturbative analysis throughout the cosmological evolution, the dilaton 
field must be bounded from above, 0(r) < 0^, with gc = e'^" sufficiently small. 



Using the state equation (4.1) and the entropy conservation (s)', valid throughout the 



cosmological evolution, we obtain 

aT = a(rj)Tc = 7oo, (i = l,...,n). (4.2) 
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The dilaton equation (iii)' can be easily integrated since p = [d — 1)P in each interval, 
yielding 



- for Ti^i <T <Ti 



for T.} < T < T, 



(4.3) 



where cf = cf,^ are constants related to the tensions (3.32): 



(4.4) 



To count the number of independent integration constants parametrizing the solution, 
we consider three independent equations, namely the iV-equation (i), the dilation equation 
(iii)' and the entropy conservation relation (s)'. These are second order in and first order 
in a, leading to three integration constants in each branch. When we glue two solutions at 
Tj, the continuity of a and across the brane, together with a(Tj_|_) = a(rj_) = have to be 
imposed. Consequently, we are left with 2x3 — 4 = 2 arbitrary integration constants. 

In the following, we display solutions characterized by a single phase transition occurring 
at cr{Tc) = 0, Tc = 0. These solutions are obtained under the assumption that Tc is reached 
only once at r = = 0. For space-time dimensions d > 2, the isotropic and homogeneous 
solutions may or may not have a non-trivial spatial curvature k. We display solutions in 
both cases. 



4.1 Bouncing cosmology with vanishing curvature, k = 

We first consider the case of vanishing spatial curvature (A; = 0), looking for a cosmological 
evolution satisfying the consistency conditions 1 to 5 listed above. As explained at the end 
of the previous subsection, the complete cosmological evolution is expected to depend on 
two arbitrary integration constants, which we may choose to be the values of the scale factor 
ttc and the dilaton field (pc at the brane. We find that the complete cosmological evolution 
in the sigma-model frame and in the conformal gauge, 

m — = m — = m — = |cr| , (4-o) 

Or ttr T 
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takes a simple form 



a 



sign(rj 
d-2 

(Pc H ; 



r7+ In 1 + 



1 



In 1 



UJac\T\ 
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ujac\T\ 

V- 

V- 



(4.6) 



where r]± = \/ d — 1 ± 1. The parameter to is proportional to the brane tension kr, which is 
responsible for the gluing: 

d-2 



Ay/d 



kr = 2^2{d-l) Tf^ 



(4.7) 



Note that kr is of order 0{e'^'' = Qc), which is larger than the naive 0{e'^'^'' = g^) expectation 



from Eq. (4.4). In the neighborhood of the brane (|/t_RacT| ^ 1), the metric is by construction 
regular while the dilaton field shows a conical singularity: 



sign(rj 



1 



d-2 

y/d^ 



1 



^^^{^Ra^rf + 0{\KRa^r\'), 









\ujacT\ 
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+ 0{{KRa,Ty). 



(4. 



The maximal value of the Ricci scalar (3.24) occurs at the bounce and it is given in terms 
of the brane tension by 

n^=^ = 0{g'^). (4.9) 



Thus both higher derivative corrections and higher genus contributions in the effective action 
remain small throughout the evolution, and can be consistently neglected, provided that the 
critical value of the string coupling, Qc, is taken to be sufficiently small. 

Far from the brane {\K,RacT\ ^ 1), the dilaton is asymptotically constant, the temperature 
drops and the scale factor tends to infinity. The whole evolution in sigma-model frame 
describes a bounce, where the scale factor and temperature are smooth. In the Einstein 
frame however, they develop conical singularities at the brane inherited from the rescaling 
with the string coupling. The Einstein frame fields are given by 



24> 



{NE,aEA/TE):= (iV,a,l/T). 



(4.10) 



At very early and late times (|a;acT| — )■ +oo), the scaling properties of the dilaton motion 
and thermal contributions to the energy density in the Einstein frame and in the Ne = 
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1 gauge are respectively l/a'^'^'^'' and The whole cosmological evolution describes 

a bounce between two asymptotically radiation dominated Universes. The fact that in 
the Einstein frame the scale factor, the temperature and the dilaton bounce with conical 
singularities follows from the localized negative contribution of the brane to the pressure, 
Pb = — e~^'^'= kr5{t). This localized pressure is not artificially added, but its origin lies in 
the stringy thermal duality properties of the system. 

To complete our analysis we also display the cosmological evolution in two dimensions: 



sign(^rj 
0" = — - — 



2 



^,-ilnfl+"^"^""" 



(4.11) 



This solution was found in the context of the Hybrid models in Ref. |16|, but it is also valid in 
the more general tachyon-free two-dimensional thermal models. Consistently, this solution 
is recovered when the space-time dimension d is formally treated as a real parameter, taking 



the limit 2 in Eqs M.6| and (4.7). 



For d > 2, figure [T] makes the comparison between (i) Classical General Relativity cou- 
pled to thermal radiation (doted lines), (ii) Classical General Relativity coupled to thermal 
radiation and non-trivial motion for the dilaton field (dashed lines), and {iii) the thermal 
superstring indicating a phase transition at Tc = between the dual effective field theories 
(sohd lines): 



(i) In the first case, the scale factor. 



'^CGR ~ 



d-2 V 



V- 



^ {ua,T + r]+ + T],) , (4.12) 
develops the well known initial singularity at infinite temperature, here at uacT = 

{ii) When the system is coupled to a non-trivial dilaton (0 oc e'^'^N/a'^~^), the cosmological 



evolution takes the form (4.6), without the absolute values in the arguments. In this 
case, cr(r) vanishes and bounces at r = 0, so that o"(r) > is always satisfied and 
the whole evolution remains within the framework of the momentum effective field 
theory. In the Einstein frame, the scale factor and temperature are monotonic and. 
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Figure 1: For k = and d > 2, comparison between: {i) Classical General Relativity coupled to radiation 
(doted lines), (ii) Classical General Relativity coupled to radiation and dilaton motion (dashed lines), and 
(Hi) the superstring picture (solid lines). Case (i) leads to the conventional initial curvature and infinite 
temperature singularities. In case (ii), the early times of the Universe are out of perturhative control. In case 
(Hi), the winding — )• momentum effective field theory phase transition at t = substitutes the previous non- 
perturbative regime with a pre-big bang cosmology at weak coupling. In Einstein frame, the whole evolution 
describes a bounce where the scale factor, the temperature and the dilaton develop conical singularities. 
Asymptotically, the two phases are radiation dominated. 



when we go backward in time, the dilaton drives the Universe into an out of control 
non-perturbative regime (0 — > +oo, qe — >■ 0, — >■ +oo when ujUcT — >■ — r^-). 

{Hi) The phase transition at r = dictated by string theory effectively substitutes the 
previous non-perturbative regime in the momentum effective field theory with a per- 
turhative pre-big bang cosmology in the dual winding effective field theory. In the 
Einstein frame, the scale factor and temperature satisfy: 



d-2-24>c 



uac\T\ 



+ 1 



u)ac\T\ 
V- 



+ 1 



a-cTc 

Qe 



(4.13) 



In two dimensions, the case (i) of Classical General Relativity coupled to thermal radia- 
tion only does not make sense. In the presence of dilaton motion, the comparison between 
cases (ii) and (iii) leads to conclusions identical to those in higher dimension. 



4.2 Bouncing cosmology with non-vanishing curvature, k 

We now turn for > 2 to the case of a homogeneous and isotropic Universe with negative 
spatial curvature k = —1. As was the case of /c = 0, the cosmological evolution is compatible 
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with the constraints 1-5 and with locahzed branes at (7{tc) = triggering the transition 
between the winding-hke and the momentum-hke field theory. The scale factor and 
dilaton 0^ at the critical temperature can be chosen as independent integration constants 
characterizing the whole solution. It is convenient to combine and (pc in terms of a 
parameter A which can be further used to define a±, 



A:= — n*i:,e^'i'^alT^ , a± 



(4.14) 



2v/^^v^rTA + v'd^ (2 + A) ' 
satisfying 1>«_|_>0 and — 1. These quantities appear explicitly in the algebraic 



expressions of the scale factor and the dilaton. In the conformal gauge (4.5), cr(r) and 0(r) 
take the following form: 



a 



sign(rj 



d-2 
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1 — a_ 
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1 — «_ 
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(4.15) 



e 2 



1 — q;_i 



For |r| ^ 1, the bounces for the metric and dilaton are respectively smooth and conical: 



sign(r) {d - 2) (2 + A) 



r' + 0(\r\ 



^(rf_l)(rf_2)(l + A) 



(4.16) 



The brane tension kc is given in terms of Oc, A and the space-time dimension, 

2^{d-l)id-2) 



VTTx. 



(4.17) 



It determines the maximal value of the Ricci scalar, which is obtained at the bounce. 



c 



0{a-') + 0{gl). 



(4.18) 



The latter is small provided the inverse scale factor and string coupling at the transition are 
chosen sufficiently small. In this case, higher derivative terms and higher loop corrections 
can be consistently neglected throughout the evolution. 

At very early and late times (|r| ^ 1), the dilaton motion vanishes while the Universe 
cools and grows to infinity. In the Einstein frame, the scale factor a^;, the temperature Te 
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ln(^^±±^^) In a 




Figure 2: For k ^ —1 and d > 2, comparison between: (i) Classical General Relativity coupled to radiation 
(doted lines), (ii) Classical Ceneral Relativity coupled to radiation and dilaton motion (dashed lines), and 
(Hi) the superstring picture (solid lines). Case (i) leads to the conventional initial curvature and infinite 
temperature singularities. In case (ii), the early times of the Universe are out of perturbative control. In case 
{Hi), the winding — >■ momentum effective field theory phase transition at t — substitutes the previous non- 
perturbative regime with a pre-big bang cosmology at weak coupling. In Einstein frame, the whole evolution 
describes a bounce where the scale factor, the temperature and the dilaton develop conical singularities. 
Asymptotically, the two phases are curvature dominated. 



and the dilaton bounce with conical behavior at the origin r = 0. In the asymptotic limits 
r —7- ±oo, the dilaton motion, the thermal radiation and the spatial curvature contribute 
to the energy density as l/a^^'^ ^\ l/a% and l/al (in the gauge Ne = !)• As a result, 
the cosmological solution describes a bounce between two asymptotically curvature dom- 
inated Universes. As a result, the cosmological solution describes a bounce between two 
asymptotically curvature dominated Universes. 



Figure [2] makes the comparison between {i) Classical General Relativity coupled to ther- 
mal radiation (doted lines), {ii) Classical General Relativity coupled to thermal radiation 
and non-trivial motion for the dilaton field (dashed lines), and {Hi) the superstring cosmol- 
ogy where two dual effective field theories are connected by a phase transition at Tc = 
(sohd lines). The conclusions are qualitatively similar to those described for = 0: 



[i) In the first case, the infinite curvature and temperature singularities occur at the initial 
time (d-2)r = ln(^^±±^^). 



,d-2 
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(4.19) 



;i - «_)'?+ V' 2 
{ii) In presence of dilaton motion (0 oc e'^'^N/a'^~^), the cosmological evolution given in 
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(4.15) with no absolute values for r behaves as follows. cr(r) vanishes and bounces at 



r = 0, which implies the evolution is described in the momentum effective field theory 
only. The dilaton, together with the Einstein frame scale factor and temperature are 
monotonic. Thus, the early times of the Universe are out of perturbative control, 
(0 — 7- +00, aE 0, Te ^ +00 when {d — 2)r — )■ lna;_,_). 

[iii) In the thermal superstring case, the phase transition at r = replaces the above non- 
perturbative regime with a perturbative pre-big bang cosmology in the dual winding 
effective field theory. In the Einstein frame, the scale factor and temperature are 



d-2 



e 2 



e 2 
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1 — Qfj 



1 — a. 



(4.20) 

ttE 



The case k = —1 gives us the possibility to consider the evolution of the Universe at 



genus-0, i.e. without the genus-1 thermal corrections. Indeed, setting A=0 in Eq. (4.17) 
the brane tension at the tree level approximation does not vanish, 

2y/{d-l)id-2) 



K 



(4.21) 



This fact indicates a non-trivial dilaton motion still exists at this order, provided d > 2. The 



explicit solution at this approximation is obtained from Eq. (4.15) in the limit A — ?■ and 
can be brought into the form. 
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where we have defined 



r 



tr 



\n{y/d^+y/d^ 



d-2 



(4.23) 



This evolution does not describe flat (Minkowski) space, as is the case of vanishing spatial 
curvature. Instead, it corresponds to a dual time-dependent curvature dominated cosmology 
which bounces at r = 0, when the spacelike brane is crossed. The fiat (Minkowski) space 
and the above k = —1 tree level cosmology are dual, as follows from the fact that they are 
both solutions at genus-0 without central charge deficit {6c = 0) of the underlying conformal 
worldsheet theory. 
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Given the fact that the notion of spatial curvature k is irrelevant in two dimensions, it 



is a non-trivial consistency check to recover the two-dimensional solution (4.11) by treating 



the space-time dimension as a continuous parameter and taking the limit d — )• 2 in (4.15). 
In particular, nc ^ i^R when ci — > 2. Since at genus-1 both thermal evolutions with k = 
and k = —1 tend to the same solution when ci — )■ 2, one expects the "dual-to-flat" classical 



cosmology (4.22) to converge to the flat (Minkowski) genus-0 phase transition in this limit. 
This is easily checked, which shows that the weak coupling limit 0c — ^ — oo and the c? — )■ 2 
limit are commuting. 



5 Conclusions 

The scope of this work was to establish a stringy mechanism able to resolve both the Hage- 
dorn instabilities of finite temperature superstring theory as well as the initial curvature 
singularity of the induced cosmology in arbitrary dimensions. The key ingredients of this 
mechanism were first isolated in the context of the two-dimensional Hybrid models whose 
right-moving sector enjoys the MSDS structure. The latter ensures boson/fermion degener- 
acy in the right-moving massive level. 

In this paper we have shown that these stringy ingredients are generic in a large class 
of = (4, 0) superstring models. Tachyon-free, thermal configurations in dimensional 

= (4, 0) models can be constructed in the presence of special "gravito-magnetic" fluxes. 
The fluxes modify the thermal vacuum by injecting into it non-trivial momentum and wind- 
ing charges, lifting the Hagedorn instabilities of the canonical thermal ensemble. The key 
property is the restoration of the thermal T-duality symmetry of the stringy thermal system, 
implying a maximal critical temperature. In all such models there are three characteristic 
regimes, each with a distinct effective field theory description: Two dual asymptotically cold 
regimes associated with the light thermal momentum and light thermal winding states, and 
the intermediate regime where additional massless thermal states appear leading to enhanced 
Euclidean gauge symmetry. 

Taking into account the genus-0 backgrounds associated with the extra massless states, we 
have shown that they source Euclidean branes localized at the critical point, which glue the 
two asymptotic momentum and winding regimes. By utilizing string calculational techniques. 
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we were able to establish that the thermal partition function can be well-approximated by 
that of massless thermal radiation up to the critical temperature. The partition function 
exhibits a conical structure as a function of the thermal modulus a, generalizing the two 
dimensional Hybrid result to any dimension. In all regimes, the equation of state is effectively 
given by 

p={d-l)P= (rf-l)n*SdT,'^e-'^l"l, 

modulo exponentially suppressed contributions. Both the energy density and pressure are 
bounded from above attaining their maximal values at the critical temperature. We show 
explicitly that the conical structure is resolved by the spacelike branes, which occur when 
cr = 0. These branes provide a localized negative pressure which is crucial in evading the 
constraints on realizing singularity- free, bouncing cosmologies, imposed by the singularity 
theorems of classical general relativity [l]. Utilizing the above ingredients, we were able to 
obtain a string effective Lorentzian action covering the three characteristic regimes simulta- 
neously. This action incorporates the spacelike brane that glues together the winding and 
the momentum regimes. 

Taking into account the localized negative pressure contribution from the branes, as well 
as the bulk thermal corrections, we have obtained non- singular analytic cosmological solu- 
tions, describing bouncing thermal and curvature dominated Universes. The bounce occurs 
at the phase transition between the momentum and winding regimes. The cosmological so- 
lutions remain perturbative throughout the evolution, provided that the value of the string 
coupling at the branes is sufficiently small. These bouncing cosmologies are the first higher 
dimensional examples, where both the Hagedorn singularities as well as the classical Big- 
Bang singularity are successfully resolved, remaining perturbative throughout the evolution. 

In this work we presented cosmological solutions associated with a single phase transition, 
with the spacelike branes occurring at a specific point in time. The brane interpretation 
leads us to believe that consistent cosmological solutions exist associated with a series of 
phase transitions, described by a distribution of branes localized in several points in time. 
At each brane the temperature reaches the maximal value. This work is currently under 
investigation |29| . 

Having in our disposal exact cosmological solutions, we can explicitly calculate the spec- 
trum of fiuctuations at early times, say at the time locations of the branes, determine their 
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propagation at latter cosmological times and compare them to the current and future obser- 
vational data. This is possible since we have analytical control on the theory describing the 
brane. This work is currently in progress 30 . 
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